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ABSTRACT 

The Kirchhoff formula for radiation from stationary 
surfaces first appeared in 1882, and it has since found many 
applications in wave propagation theory. In 1930, Morgans 
extended the formula to apply to surfaces moving at speeds 
below the wave propagation speed; we refer to Morgans’ 
formula as the subsonic formulation. A modem derivation of 
Morgans’ result was published by Farassat and Myers in 
1988, and it has now been used extensively in acoustics, 
particularly for high speed helicopter rotor noise prediction. 
Under some common conditions in this application, however, 
the appropriate Kirchhoff surface must be chosen such that 
portions of it travel at supersonic speed. The available 
Kirchhoff formula for moving surfaces is not suitable for this 
situation. In the current paper we derive the Kirchhoff 
formula applicable to a supersonically moving surface using 
some results from generalized function theory. The new 
formula requires knowledge of the same surface data as in 
the subsonic case. Complications that arise from apparent 
singularities in the new formulation are discussed briefly in 
the paper. 


1. INTRODUCTION 

Since its first publication in 1882, the Kirchhoff formula 
has played a fundamental role in the study of optical, acoustic 
and electromagnetic wave propagation. 1 " 2 The formula 
provides a representation of solutions to the homogeneous 
wave equation in regions interior or exterior to a closed 
surface in terms of data specified on the surface itself. Here 
we will call this surface the Kirchhoff surface. The original 
formula applies only to a stationary Kirchhoff surface. In 
many circumstances, however, it is convenient to have an 
analogous formula applicable to a moving surface, and such 
a result was derived by Morgans, 3 who extended the original 
formula of Kirchhoff to include surfaces moving at speeds 
below the wave speed in the propagation medium; we refer 
to Morgans’ result as the subsonic Kirchhoff formula. A 
modem derivation this formula based on generalized function 
theory and in a form suitable for application was given by the 
authors in 198 8. 4 

In addition to the fact that representation formulas of the 
Kirchhoff type are of interest from a fundamental analytical 
standpoint, they have in recent years also assumed importance 
as computational tools for use in conjunction with CFD 
calculations. One illustration of the utility of the subsonic 
Kirchhoff formula is its successful application by Lyrintzis 5 
to helicopter noise prediction. Here unsteady aerodynamic 
calculations were performed in the near field of a helicopter 
rotor in a reference frame fixed to the rotating blades (see 


fig. 1). These then provided data on a blade-fixed Kirchhoff 
surface like that indicated in fig. 1 for subsequent application 
of the subsonic Kirchhoff formula to determine the noise 
radiated by the blade. This particular application, however, is 
limited by the fact that under some commonly occurring 
conditions the shock system associated with a high speed 
helicopter rotor can extend well beyond the tip region, a 
phenomenon that has been called delocalization. 6 To apply 
a (linear) Kirchhoff formula in this case requires that the 
entire shock system be included inside the Kirchhoff surface. 
Because part of the surface thus travels at supersonic speed, 
the available Kirchhoff formula for moving surfaces is not 
suitable. 

The primary difficulty with the subsonic formula is the 
appearance of the Doppler factor 1-M,. in the denominator of 
the integrands, where M r is the component of the surface 
Mach number in the radiation direction. For supersonic 
surfaces, there exist directions in which M r =l, and thus 
singularities appear in the subsonic Kirchhoff formula. Here 
we derive a new Kirchhoff formula that does not contain the 
Doppler singularity. The new formula is actually valid for all 
surface speeds, but, because it is somewhat complex to code 
for computer applications, its use is not recommended for 
Kirchhoff surfaces, or portions or Kirchhoff surfaces, moving 
subsonically. In this paper, therefore, we refer to the new 
formula as the supersonic Kirchhoff formula. The analysis is 
carried out with a view toward practical numerical 
implementation. Thus, the new formula is derived for an 
open surface panel considered as a portion of the closed 
Kirchhoff surface. This allows its use to be restricted to just 
those parts of the Kirchhoff surface that are actually moving 
at supersonic speed. 

In the next section, the inhomogeneous source terms of 
the wave equation leading to the Kirchhoff formula are 
derived. In the following section, the solution of this wave 
equation is obtained which is valid for supersonically moving 
surfaces. This solution is the desired Kirchhoff formula. In 
the subsequent section the singularities associated with the 
new formula are briefly discussed. The main complications 
in the supersonic formulation are the existence of multiple 
emission times and the appearance of small or vanishing 
denominators in the algebraic expressions. We indicate how 
these complications can be overcome in applications. 

2. THE INHOMOGENEOUS WAVE EQUATION 

Our approach to deriving the supersonic Kirchhoff 
equation is similar to that of reference 4. The primary 
reference for the mathematics used in the derivation is a 
recent publication by the first author. 7 Let us assume that the 
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moving Kirchhoff surface is defined byf(x,t)=0 with f>0 
in the exterior of the surface and f<0 in the interior. This 
function is selected in such a way that |Vf|=l over the 
surface, which can always be done. 7 We are interested in 
radiation into the region exterior to the surface f=0. Let 
<J>(x,t) satisfy the homogeneous equation in this region: 

-^$>=0 (1) 

c 2 at 2 


in which c is the speed of sound. To find the inhomogeneous 
source terms of the wave equation leading to the Kirchhoff 
formula, we extend <t> to the entire three-dimensional space 
as follows: 


*(x,t) 



f>0 

f<0 


( 2 ) 


We now find I if O where D 2 stands for the wave operator 
with generalized derivatives. 

Using the rules of generalized differentiation, 7 ' 10 we 
obtain 

M® 6 (f) (3-a) 
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us assume that a panel on f=0 is defined by the equation of 
its edge curve, f=0. We assume that f>0 on the panel. We 

can always define f in such a way that Vf=v where v is 
the local unit inward geodesic normal at the edge of the 
panel. The geodesic normal is tangent to the surface f=0 and 

is normal to the edge of the panel f =f =0 (see Fig. 2-b). For 
the open surface at this panel, the source terms of the wave 
equation in eq. (5) must be modified by the Heaviside 

function H(f) as follows: 

ff$ = _[^!l^. + $ ii |H(f)6(f) 

1 c 3t ) (6) 
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Now we carry out the time differentiation and divergence in 
the last two terms of eq. (6). We have 


[M <D H(f) 6(f)] = i (M„ <5)H(f)6(f) 
c dt c ot 

-M n M v 4> 8(f) 6(f) -M 2 fcH(f)5'(f) 


where a tilde under a function stands for the restriction of the 
function to the surface f=0. 7 The symbol M v =v v /c is the 

local Mach number of the edge f = f = 0 in the direction of the 
geodesic normal. Similarly, the last term in eq. (6) can be 
written as 


In eqs. (3) the bar over a derivative stands for generalized 
differentiation, and M n =v n /c is the local normal Mach 
number of the surface f=0. The Dirac delta function with 

support on f=0 is denoted 5(f). The generalized Laplacian of 
is similarly found as follows: 

V 6 = V<t> + <£nS(f) 

V 2 6=V 2 «i>+<I> n 6(f)+V-[«l>a5(f)] (4_ b ) 


Here n = Vf is the local unit outward normal and 
= n -V^> on f =0. From eqs. (3) and (4), and the fact that 

=0, we find 


n 2 ® 
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^ c dt 


6(f) -I-f [M n O 6(f)] -V-[*fi5(f)] 
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(5) 


The last two terms of this equation need further 
manipulations to make them suitable for the derivation of the 
supersonic Kirchhoff formula. We perform these 
manipulations below. 

In applications, the Kirchhoff surface is first divided into 
panels (Fig. 2-a). Some of these panels travel at supersonic 
speed, so that the subsonic Kirchhoff formula is not suitable. 
We therefore, derive the new formula for a single panel. Let 


V-[^nH(f)6(f)] = -2H f 0>H(f)6(f) + $H(f)6 / (f) (8) 


where H f is the local mean curvature of the surface f=0. 
After substituting eqs. (7) and (8) into eq. (6) and collecting 
terms, we get 


rf<S> = - 0) + — — +-~(M 0)-2H f 0 

c at car 1 ' f 

- (1 -M 2 ) $H(f) 6 '(f) +M n M v $ 6(f) 5(f) 


H(f)6(f) 
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This is the desired form of the inhomogeneous source terms 
of the wave equation. In the next section, we give the 
solution of eq. (9) which is the supersonic Kirchhoff formula. 


3. THE SUPERSONIC KIRCHHOFF FORMULA 


Equation (9) is now written as 

□*$ =q,(x,x)H(f)6(f) + q 2 (x,T)H(f)6'(f) 
+ q 3 (x,T)6(f)5(f) 


Note that q 2 is restricted to f=0 and hence is written as q 2 . 
Next we write ® as the sum of three functions 
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o =3>J +$ 2 + C >3 


where these functions are the solutions of the wave equations 

S J 3. 1 =q 1 (x,t)H(f)5(f) (12 _ a) 

i=f« 2 -g 2 G?.t)Htf)6 / (Q (12-b) 

Cf<t»3=q 3 (x,t)8(f)5(f) (12_C) 


Pi _* n-M„r VF 

— =N-V , N = —2- = 7~t 

aN A JVF| 

cos6 / =N-N , ) $'--1^1--^- 


A 2 = 1 +M 2 -2M„cos0 


cos 0 = v * r 


(16-a,b) 
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The solutions of these equations are given fully in reference 
7 (see solutions of eqs. 4.23-b, e, f in Ref. 7). Here we 
utilize these results directly and refer readers to the original 


The solution of eq. (12-a) is given by eq. (4.34) of 
reference 7 as follows; 

4n 3>j(x,t) = f ^dE (13) 

F=0 lA 
F>0 

where Qj = [q 1 (y,'r)] ret =q 1 (y,t-r/c) , r=|x-y|, 
F(y;x,t)=[f(y,T)] ret , F(y;x,t) =[f(y,T)] M and 

A 2 = 1 +M^ -2M n cos0 ( 14 > 


In eq. (15) we have also denoted the mean curvature on the 
S-surface by H F , and dL is the element of the length of the 
edge of the Z-surface. Below we will see why we must 
retain the restriction sign (~) in Q 2 in the integrand of the 

surface integral. 

The solution of eq. (12-c) is given by eq. (4.49) of 
reference 7 as follows: 


4it* 3 (X,t)= f -2idL 

rA, 


F=0 *A> 
F=0 


where Q 3 = [q 3 (y,t)] ret and 


A o =|VFxVF|-AAsin0 / 


Here M n is local normal Mach number on the panel and 
cos 0 — n * f , where f=(x- y )/r is the unit radiation vector. In 
eq. (13) d£ is the element of surface area of the surface 
F = 0,F>0. This is the surface formed by the curves of 
intersection of the collapsing sphere r =c(t“i), (x,t) fixed, 
with the panel in motion as the source time t increases from 
to t. Figure 3 illustrates the construction of the S-surface for 
a panel. 

The solution of eq. (12-b) is given by eq. (4.42) of 
reference 7 in the form 

4ti 0 ? (x,t) = / (— ) + d£ 

2 p { 0 [ A3NV rA 2 j 

-/ ?- 2 --^dL (15> 

J rA 2 
p=n m 


where Q 2 = [q 2 (y , t )] ret and we have introduced the following 
symbols: 


After combining the eqs. (13), (15) and (17) we get 
. f 1 L 2H f n N-VA„ N-VqJ 


f i^r Q dE + f J_L Aoob^ l d (19) 

L r 2 A 2 2 L 3 A 2 J 


In eq. (19) we have explicitly performed the derivative d/5N 
in eq. (15) and separated the near and far field surface 
integrals. We note that 

N-VQ 2 =[N-VqJ --[N-VrqJ ( 20 ) 


where the symbol q 2 is used for dq 2 (y 9 x)/dx . From eq. 
(16-b), we have 

t 1 

N^[n-M n (ncos0 +t 1 sin0)j = -^(l -M n cos0)n 
M sin0 _ (21) 


in which 1 1 is the unit vector along the projection of r on 
the local tangent plane of fr=0. Since dq 2 /3n=0 (this a 

property of the restriction 7 ), we get 
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tf*Vq 2 


M a sin8 3q 2 (y >x) 
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where d/do x is the directional derivative in the direction of 
t j (keeping x fixed). We also have 


N *Vr = 


M„ - cos 0 

n 

A 


(23) 


Kirchhoff formula from eq. (25) when the surface is 
stationary. We assume a closed surface described by 
f(y)=0. In this case there is no line integral and we have 


F(y;x,t) 

= f(y)=0 


(26-a) 

M„=0 , 

A =1 


(26-b,c) 

H P =H f . 


-2H f <5> 

(26-d,e) 

= $ - 

q 3 =0 , 

dE=dS 

(26-f,g,h) 


which, when substituted in eq. (20), gives 


N-VQ 2 * 


M sinO dQ 2 cos0-M 


do. 


cA 


Q 2 (24) 


In eq. (24), Q 2 stands for (dq 2 /dx) rc< and the directional 
derivative again is calculated keeping x fixed. Finally, after 
using eqs. (23) and (24) in eq. (19), we get 


4tt <i> (x,t) = f — 
J r A 


F-0 

P>0 
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2H P + N-VA'I 

rA 
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A + A 2 J 


M n sin0 0Q 2 cos0~M a 
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Equation (25) is the supersonic Kirchhoff formula and is the 
main result of this paper. We will discuss the above result 
further below. 

First, however, we remark on the significance of utilizing 
the restriction of the function q 2 to the surface f=0 in the 
derivation of eq. (25). The key property of the restriction is 

that dq 2 /dn = 0, which allows the derivative in the direction N 

in eq. (15-b) to be evaluated very simply (see eq. (22)). If 
we had not introduced the restriction, the ultimate result 
would still be eq. (25), but it would contain a number of 
extra terms arising from unnecessary differentiation of q 2 in 
the direction n. The extra terms, of course, cancel one 
another, but is difficult to recognize the cancellation in the 
form in which the terms arise in eq. (25). Using q 2 leads 

directly to the above result which the authors believe is the 
simplest possible expression of the supersonic Kirchhoff 
formula. 


where dS is the element of the surface area of f=0. 
Substitution of these results in eq. (25), yields 

4 7 i <I>(x,t) « j — [c^cosO <i> - O n ]^dS + J flcosQ ^ 

f-0 r f=0 r 

(27) 


This is the classical Kirchhoff formula.’ 2 Thus this aspect of 
consistency of our main result has been established. It is 

evident that eq. (25) depends on and <i> on the 

Kirchhoff surface f=0. This is, of course, expected. Note 
that dQ 2 /do x involves d®ldo p which is obtained from 
knowledge of <D on f=0. 

We will now discussthe problem of singularities of eq. 
(25). The quantity A appears in every term of the 
denominator and produces a singularity when A=0. The 
singularity ^=0 appears in the line integral of eq. (25) if 
the conditions discussed in reference 1 1 are satisfied. When 
A=0, the collapsing sphere is tangent to the Kirchhoff surface 
f-0 at a point where M n =l. The E surface in the vicinity of 
this point is very complicated with a possible mean curvature 
singularity. This case requires further analysis which we will 
not pursue at present. For now we will indicate a practical 
method to get around this problem. Assume that the 
Kirchhoff surface is chosen to have a shape like that of a 
biconvex airfoil such that there are no points on the surface 
at which M n =l. Then we only have to be concerned with the 
singularities of the line integral in eq. (25) on the leading and 
trailing edges of the Kirchhoff surface. However, the 
integrand of the line integral of eq. (25) is precisely eq. (17- 
c) of reference 1 1 with p'c 2 replaced by <t>. An analysis 
identical to that of reference 1 1 shows that the singularity of 
the line integral in eq. (25) is integrable so long as is 
continuous in the vicinity of the singular point. But this 
condition is always satisfied in practice so that there will be 
no numerical problems in using eq. (25) for the proposed 
Kirchhoff surface. 


One important fact that we point out here is that eq. (25) 
is valid for a Kirchhoff surface in arbitrary motion, including 
motion at subsonic speed. However, we do not recommend 
its use for subsonic speeds because the formula presented by 
the authors in reference 4 is much more efficiently applied 
for computation. 

4. ANALYSIS OF THE MAIN RESULT 

We first check to see if we recover the classical 


5. CONCLUDING REMARKS 

In this paper we have derived a Kirchhoff formula that 
describes solutions of the homogeneous wave equation 
exterior to a surface in arbitrary motion. The formula is 
specifically designed for practical computation of radiation 
from surfaces in supersonic motion. Because surfaces 
moving at subsonic speeds are more efficiently treated using 
an earlier result, the new formula is derived for an open 
surface; it is recommended that it be used only on surface 
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panels that are actually moving supersonically. The new 
formula is expressed in a relatively simple form in terms of 
surface data and of easily calculated geometric and kinematic 
properties of the moving Kirchhoff surface. It is complicated 
somewhat by the existence of singularities that occur under 
certain conditions in supersonic motion, but these singularities 
are shown to present no difficulties in numerical 
implementation of the formula. 
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Fig. 2 Open panel on Kirchhoff surface. 



Fig. 3 Construction of E-surface for panel. 
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